Macroscopic Hong-Ou-Mandel interference 
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We report on a Hong-Ou-Mandel (HOM) interference experiment for quantum states with photon 
numbers per mode as large as 10^ generated via high-gain parametric down conversion (PDC). The 
standard technique of coincidence counting leads in this case to a dip with a very low visibility. By 
measuring, instead of coincidence rate, the variance of the photon-number difference, we observe 
a peak with 99.99% visibility. From the shape of the peak, one can infer information about the 
temporal mode structure of the PDC radiation, including the degree of frequency /time entanglement. 

PACS numbers: 42.50.-p, 42.50.Lc, 42.50.Ar, 42.65.-k 



One of the most fascinating manifestations of quantum 
interference is the Hong-Ou-Mandel (HOM) dip, an ef- 
fect that started a whole new direction in quantum optics 
First reported in 1987 [H [H, it consists of the interfer- 
ence of two single-photon wavepackets on a 50% beam- 
splitter (BS). When indistinguishable photons simultane- 
ously arrive at different input ports of a BS, they always 
depart from the same port, due to the destructive inter- 
ference of the probability amplitudes to be both reflected 
or both transmitted The effect is usually observed 
as a 'dip' in coincidences, its width determined by the 
coherence time of the photons and not affected by the 
group- velocity dispersion Q . HOM interference is widely 
used for identity tests of single-photon states @ , includ- 
ing ones generated by a quantum dot and a laser 
conditionally prepared via parametric down-conversion 
(PDC) and a highly attenuated coherent state 0], pro- 
duced via two consequent emissions of a single atom Q 
or molecule Q or by two different molecules [l^l- Other 
important applications are the implementation of a pho- 
ton Fock-state filter [llj and quantum-optical coherence 
tomography 12]. It is worth mentioning that if, as in 



most cases, the photons are generated via PDC, the level 
of coincidence rate in the 'dip' is not zero but is given by 
the accidental coincidences rate [l^l ■ For this reason, the 
visibility of the dip reduces down to 33% at high-gain 



PDC jl4l | and becomes vanishingly low if many modes 



are involved. 

In this paper we expand the boundaries of the HOM 
effect applications and observe destructive HOM interfer- 
ence for bright twin beams (containing, on the average, 
up to 8 • 10^ photons per mode) generated via high-gain 
PDC in a traveling-wave optical parametrical amplifier. 
Namely, instead of two single-photon states we sent to 
the beamsplitter two macroscopic states containing ex- 
actly equal number of quanta. We suggest an alterna- 
tive approach of observing the HOM effect, based on the 
measurement of the normalized variance of the difference 
signal at the output ports of the beamsplitter. This tech- 
nique is robust against the multi-mode detection fl5l | and 



allows one to observe the interference with a high visibil- 
ity. 

The experimental setup (Fig [1]) contained three main 
parts. In the state generation part, horizontally polar- 
ized third harmonic of a pulsed Nd:YAG laser with the 
wavelength 354.7 nm, pulse duration 18 ps, and repeti- 
tion rate 1 kHz was used as a pump. The intensity was 
changed by rotating a half-wave plate (Ap/2) in front 
of a Glan-Lascr Polarizer (CP). The bright two-mode 
polarization squeezed vacuum state was generated in a 
frequency-degenerate traveling-wave optical parametric 
amplifier (OPA) based on two 5-mm width crystals BBO 
cut for type-II phase matching. The effect of the spatial 
walk-off was reduced by placing two crystals with the 
optical axis in the horizontal plane tilted oppositely with 
respect to the pump. The OPA was tuned to generate 
twin beams in the collinear direction. A high parametric 
gain was achieved by reducing the pump beam diame- 
ter to 180 with the help of a telescope consisting of 
a convex lens {F = 50 cm) followed by a concave lens 
{F = —7.5 cm) at a distance of 42.5 cm. Due to the 
group velocity difference, the ordinarily polarized pulse 
of the signal radiation was always delayed relative to the 
extraordinarily polarized idler pulse. After the crystals 
the pump was eliminated by two dichroic mirrors DMi_2 
with high reflection at 709.3 nm and high transmission at 
354.7 nm. The rest of the pump was absorbed by the red 
glass filter RG-630 with AR coating in the spectral range 
of the PDC radiation. The delay between the signal and 
idler beams was changed in the delay line with the polar- 
ization beamsplitter (PBSi) at the input by scanning the 
position of the mirror Mi with respect to the mirror M2. 
The linear polarization of the signal and idler beams was 
rotated by 90 degrees after double passing through the 
quarter wave plates (As/4) oriented at 45 degrees to the 
PBSi. Therefore both beams left the interferometer from 
the same output and were directed to the registration 
part of the setup. A narrow angular spectrum (with the 
width 4 mrad) was restricted by a 3 mm iris aperture (A) 
placed in the focal plane of a 75 cm collecting lens (Lg). 
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The HOM interference was observed at the output of the 
polarization beamspUtter (PBS2) placed after a halfwave 
plate (As/2) which rotated the polarization basis of the 
PBS2 by 45 degrees. In this case, one can show that 
the trajectories of two multi-photon states arriving si- 
multaneously to the polarization beamsplitter became in- 
distinguishable and the interference occurred. After the 
beamsplitter the radiation was focused by two lenses (L) 
on the analog detectors (Di^2) based on pin-diodes jl6l |. 
The electronic pulses from the detectors had a fixed du- 
ration (8^s) and their amplitudes were proportional to 
the numbers of photons at the input. Then they were 
time- integrated with the help of a 10-Msample/s 14-bit 
analog-to-digital converter card. The card was triggered 
by the synchronization pulses of the laser. The measured 
signals per pulse 81^2 at the output of D 1^2 were used for 
the calculation of the normalized second-order intensity 
correlation function g^^^ = and the normalized 

variance of the difference signal ttt-t = ^t^-tttS^- The 

o {Si} + {S2) 

angle brackets denote the averaging over an ensemble of 
pulses. In the experiment, the variance of the difference 
signal, the second-order intensity correlation function, 
and the mean values of the signals were calculated for 
each position of the mirror Mi by averaging the data 
over 30000 pulses. 

At the beginning, the brightness of the generated state 
was evaluated from the parametric gain (G) measure- 
ment. The intensity of the PDC radiation was measured 
as a function of the pump power (see bottom-left part 
of FigllJ. The measurement was done after narrowband 
spectral selection (0.2 nm around the wavelength of 709.3 
nm) by means of a spectrometer HORIBA Jobin Yvon 
Micro HR (not shown) and angular selection (4 mrad 
around the coUincar direction) by means of an aperture. 
The red circles represent the experimental results and 
the blue line is the fit plotted according to the depen- 
dence N = sinh^ G. The curve clearly demonstrates that 
the process was strongly nonlinear. The gain reached the 
value of G = 7.5 at 55 mW pump power, which is equiv- 
alent to Nmodo ~ 8 • 10^ photons per mode. Fig. [2] (a) 
shows the interference dip in the second-order intensity 
correlation function g'^^ plotted versus the delay time. 
According to the definition of the interference visibility 
V = "°^7""" , the obtained value of ^ = 0.022 was very 
small. The value of the g*-^^ at the edges was used for 
the estimation of the number of detected modes. In this 
case, the relative delay was large enough and there was 
no interference at the beamsplitter. Hence, the initial 
correlation between the signal and idler beams given by 
the second-order intensity correlation function for a sin- 
gle mode g*-^' = 2-1- jy ^ ^ was conserved. Assuming that 
due to the multi-mode detection the initial normalized 
second-order intensity correlation function g^'^^ is trans- 
formed into gm^ = 1 + ^ , the number of detected 
modes was found to be m = 10. 




FIG. 1. (Color online) Experimental layout for the macro- 
scopic HOM interference: Nd:YAG (3a;), third harmonic of 
the Nd:YAG laser; Ap/2, half-wave plate for 354.7 nm; GP, 
Glan prism; M, mirror; BBO, two 5 mm type-II beta bar- 
ium borate crystals with the optical axes tilted oppositely 
with respect to the pump; DM1,2, dichroic mirrors with high 
transmission for the pump and 99.5% reflection for the down- 
converted radiation; KG, red glass filter RG630 with AR coat- 
ing at 709.3 nm; Ls, lens with the focal length 0.75 m; PBSi,2, 
polarization beamsplitter; As/4, quarter-wave plate for 709.3 
nm; Mi, 2, high-reflection mirrors at 709.3 nm; A, 3 mm aper- 
ture placed in the focal plane of the lens Ls; As/2, half- wave 
plate for 709.3 nm; L, lens; Di,2, detectors. Bottom left: The 
dependence of the PDC signal on the pump power. 



The normalized variance of the difference signal is dis- 
played as a function of the time delay in Fig. [5] (b). Sim- 
ilar to the two-photon HOM interference, multi-photon 
signal and idler states arriving simultaneously to the 
beamsplitter were mostly directed as a whole either to 
one port or to the other. Therefore, the output states 
of the PBS2 were anticorrelated leading to the observa- 
tion of a dip in the correlation function or to a peak in 
the variance of the difference signal. Indeed, a narrow 
central peak was observed in the experiment. The width 
of this peak is determined by the coherence time of the 
PDC radiation. The broader pedestal peak has a classi- 
cal nature. It results from the pulsed shape of the signal 
and idler radiation. One should mention here that due 
to the strong nonlincarity of the parametric process the 
pulse duration of the PDC radiation was much shorter 
than the pulse duration of the pump. As it was dis- 
cussed above, the signal and idler pulses arriving to the 
beamsplitter in sequence were split independently and 
their excess fluctuations were eliminated as a result of 
subtraction. Therefore the value of the normalized vari- 
ance of the difference signal at the edges corresponds to 
the shot noise level. As a result, the peak manifests a 
visibility V — 0.999 exceeding all previous records of 
HOM interference, including the one achieved recently 
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FIG. 2. (Color online) HOM interference: dependencies of the 
second-order intensity correlation function (a) and normalized 
variance of the difference signal (b) on the delay between sig- 
nal and idler beams. Red circles represent the experimental 
results and the blue line is the theoretical dependence. 



with two oppositely chirped classical optical pulses in an 
experiment on sum frequency generation [l9j . 

The theoretical curve shown by the blue line is plotted 
according to the approach discussed in detail in [17| . Here 
we use the Hcisenberg picture of a quantum state evolu- 
tion within the undepleted pump approximation. In this 
case, the Hamiltonian describing coUinear typc-II OPA 
with the signal and idler frequencies (±17) tuned around 
the degenerate frequency ^ is 



H oc T{n)[a\{n)al{-n)] + H.c, 



(1) 



where T{n) is the spectral parametric gain coefficient and 
a\ , a2 are the photon creation operators in the signal and 
idler beams. Considering the coUinear geometry of the 
experiment the calculations are done within the plane- 
wave approximation. The state at the output of the OPA 
is described using the Bogoliubov transformations [3| 



a2iil) = U{n, i)ar'(f^) + Vi^^ t)al'"''' (^Q), (2) 



where al'"',a'^° 



,02 



are the photon annihilation 



and creation operators at the vacuum inputs of the OPA. 
The gain functions U{n,t) and V{i},t) are given by 



U{n, t) = coshr(f], t) + ^ A|^sinhr(l], i), 



(3) 



with 



(4) 



where Ic is the total length of the crystals and G{t) de- 
pends on the the pump field envelope ^o(0 ^(0 = 
aAo{t)lc- Here ct is a coefficient depending on the prop- 
erties of the crystal. As long as the frequency spec- 
trum of type-H PDC is narrow (1.3 nm), the phase mis- 
match A is defined only by the temporal walk-off be- 
tween the signal and idler waves at the degenerate fre- 
quency A{n) ^ [{dki/duj)^^^j2 ~ {dk2/duj)^^^^/2]fl. 
Assuming that the signal and idler states arrive at the 
beamsplitter with the relative delay r, the beamsplitter 
transformations take the form 



71 ' ' 



-ai{t - r) + 02(0 
71 



ar*(i) 

(5) 

Notice that if the pulse duration is Tp the number of 
photons per pulse at each output of the beamsplitter is 

A^i.2 = /2»r2^W»r2 W^*- According to the defi- 
nition of the normalized variance of the photon-number 

difference ^ ^ ^^^^ ^^^^^ 

brackets denoting the averaging over the vacuum state, 
we obtain the expression for the numerical evaluation, 



(iv+) /o°°df7|y(o,i7)|2 
dr!|i/(o, n)\''{\v{T, rj)|2 + rc([c/*(o, n)]''e''^^^)}. 

(6) 

Taking into account the finite quantum efficiency -q of 
the detection, the actual normalized variance of the dif- 
ference signal is 



A^^(r) 



= 1 



, AiV^(r) 



(7) 



In accordance with the experiment, numerical calcula- 
tions were done for two type-II BBO crystals with the to- 
tal length of Ic = 10 mm, assuming that the pump pulse 
has a Gaussian envelope with the duration Tp = 18 ps, 
and a parametric gain of G = 7.5. The obtained value 
of the fitting parameter is 77 = 0.03. The result can be 
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explained by the losses in the optical channel and by the 
mismatch of the signal and idler beams. This mismatch is 
due to the fact that although signal and idler beams have 
different angular spectra because of the anisotropy, they 
are restricted by the same aperture A. Therefore, we find 
that the theoretical estimations are in a good agreement 
with the experiment. It is also worth mentioning that 
the central narrow peak is twice higher then the broad 
one and its amplitude is twice the mean photon number 
per mode 2 • N. Since the background is given by a unity, 
we conclude that the brighter the state the higher visi- 
bility can be obtained. The ratio between the widths of 
the broad pedestal and the central narrow peak gives the 
effective number of the longitudinal modes, miong = 8. 
It is important to mention here that for pure states this 
ratio is an estimate for the degree of entanglement in the 
time/frequency domain pH . 

Figure [3] shows the full width at half maximum 
(FWHM) of the central peak as a function of the para- 
metric gain. A 20% narrowing of the peak was ob- 
served by changing the parametric gain in the range of 
5.5 < G < 7.5. The red circles are the experimental 
data, in perfect agreement with the theoretical predic- 
tion (the blue line) . As it was shown in , at high gain 
most of the photons are generated at the output part of 
the crystal, which leads to the broadening of the spec- 
trum and therefore, to the reduction of the coherence 
time of the PDC radiation. This makes the experiment 
more sensitive to the balancing of the arrival times of two 
multi-photon states and therefore significantly increases 
the accuracy of the time measurement. 




5.5 6.0 6.5 7.0 

Parametric gain 



7.5 



FIG. 3. (Color online) Full width at half maximum (FWHM) 
of the HOM peak as a function of the parametric gain. Red 
circles represent the experimental results and the blue curve 
is plotted according to the theoretical model. 

To summarize, we have observed the Hong-Ou-Mandel 
interference of macroscopic states containing on the av- 
erage 8 • 10^ photons per mode. By measuring the nor- 
malized variance of the difference signal we were able to 
achieve the visibility of 0.999. This is an evident advan- 



tage over the correlation function measurement where 
the visibility is strongly reduced in the presence of high 
gain and many modes. The experimental results provide 
the information on the temporal mode structure of the 
interfering fields. Moreover, we believe that this mea- 
surement can be used for the estimation of the degree of 
frequency/time entanglement for pure states. It is shown 
that at high gain PDC the width of the interference peak 
reduces with the parametric gain leading to an increase 
in the time resolution. Therefore, we can conclude that 
the HOM effect for the bright states generated via high- 
gain PDC has the same properties as the HOM effect for 
single photons but a better resolution. 
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